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We show that repulsive local Coulomb interaction alone can drive valence-skipping charge dis- 
proportionation in the degenerate d-band, resulting in effective negative-U. This effect is shown to 
originate from anisotropic orbital-multipole scattering, and occurs only for d 1 , d 4 , d 6 and d 9 fillings 
(and their immediate surroundings). Explicit boundaries for valence-skipping are derived and the 
paramagnetic phase diagram for d 4 and d 6 is calculated. We also establish that the valence-skipping 
metal is very different, in terms of its local valence distribution, compared to the atomic-like Hund's 
metal. These findings explains why transition metal compounds with the aforementioned d-band 
fillings are more prone to valence-skipping charge order and anomalous superconductivity. 
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While going up in an elevator, have you ever caught 
yourself staring at the numbers flying by? The eleventh, 
the twelfth, and then all of a sudden the fourteenth floor. 
Valence-skipping elements are just like elevators, except 
they avoid certain valence states rather than the thir- 
teenth floor. The most prominent valence-skippers are 
the post-transition metals, Tl, Bi, Sb, etc., who display 
missing valences in many of their compounds pQ. In 
general valence-skipping is driven by a negative effective 
Coulomb repulsion U e g, but the mechanism causing this 
is debated. Anderson [2] showed that static lattice re- 
laxation can drive U e g < 0. However Varma pQ noted 
that even free atoms have reduced U c e in the unfavor- 
able valences, and proposed an intra-atomic electronic 
mechanism. But more recently the electronic route has 
been discredited for these elements, in favor of the lattice 
relaxation mechanism [3|. 

Apart from the post-transition metals, valence- 
skipping has also been observed in transition metals p] 
with a d n — > d n ~ 1 +d n+1 type of charge disproportion- 
ation. Experimentally the most evident examples are 
the iron compounds (La,Ca)Fe03 [H |S], (La,Sr)Fe03 
[51 ©, and Sr 3 Fe 2 07 [7., where Mossbauer spectroscopy 
has established valence-skipping Fc 4+ — » Fe 3+ +Fe 5+ , 
(d 4 — > d 3 +d 5 ) charge-order, even in the absence of lattice 
relaxation [?J [7J. Theoretically Katayama-Yoshida and 
Zunger [5] showed that effective monopole screening of 
intra-atomic interactions indeed can give rise to U e s < 
in transition metal impurities. This idea has been used 
to explain the charge-order in YN1O3 [9] within a two- 
band e g model. But the complete d-band still deserves 
more attention. 

The fact that some authors even refer to valence- 
skipping as "mysterious" [TU1 HJ, show the need for 
better understanding of the underlying mechanism be- 
hind this phenomenon and its systematics. The result- 
ing negative-U model however, has been studied exten- 
sively, and shown to drive both charge-order and anoma- 
lous superconductivity 12J. So unveiling the mysteries 



of valence-skipping could pave the way for more exotic 
physics. 

In this Letter we show that valence-skipping in the de- 
generate d-band is driven by the higher orbital-multipole 
part of the repulsive intra-atomic Coulomb interaction. 
This effect is found to be limited to the particular fillings 
d 1 , d 4 , d 6 , and d 9 , and their immediate surroundings. 
From a multiplet analysis we derive explicit bounds for 
valence-skipping, and finally the emerging anomalous va- 
lence fluctuations in the paramagnetic metal are studied. 

Let us begin by constructing a minimal model for the 
correlated d-band. We assume that the Coulomb inter- 
action is local and rotationally invariant, a good first ap- 
proximation for transition metals |13j . Under this as- 
sumption, the interaction is exactly given by the Slater- 
Condon angular-momentum expansion, and the Slater- 
integrals F (0) , F (2) and F (4) [14J. For the electron hop- 
ping we use a degenerate semi-circular density of states, 
and take the half-bandwidth as our unit of energy. 

In general the local interaction describes electron-pair 
scattering between local two-particle states, and rota- 
tional invariance ensure that these processes conserve lo- 
cal total orbital momentum L and total spin S. As we 
are going to see, anisotropic orbital-multipole scattering 
(i.e. for L > 0) has an intrinsic connection to valence- 
skipping. To make this clear we now seek to isolate this 
contribution to the interaction. 

Within the Slater-Condon interaction, the term 
is a density-density interaction with isotropic scattering, 
while the F^ and F^' terms have different scattering 
strengths for all L and S. Interestingly their orbital- 
multipole anisotropies cancel out when F( 4 ) /Ja 2 ) = 9/5. 
This corresponds to a Laporte-Platt degenerate point of 
the Slater-Condon interaction with large accidental de- 
generacies of multiplets [15] . 

Spurred by this observation we investigated the Slater- 
Condon interaction in detail [16] , and found that in this 
point it simplifies to the rotationally invariant Kanamori 
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interaction [17] . having the compact form 

H = (U- 3J)N(N -l)/2 + J{Q 2 -S 2 ), (1) 

where N is the total number operator, and S 2 and Q 2 
are the total spin and quasi-spin operators [14] . 

The Kanamori interaction and the roles of its coupling 
parameters, the Hubbard U, and Hund's rule J have al- 
ready been studied extensively [TB]- This makes the re- 
duction of the Slater-Condon interaction (at the Laportc- 
Platt degenerate point F (4) /F^ = 9/5) very interesting. 
At this point U and J alone can be used to determine 
F (k \ 

But let us first establish our claim that anisotropic 
orbital-multipole scattering is indeed missing in Eq. ([lj. 
We need not to worry about the density-density interac- 
tion giving isotropic scatteringflike the F^ term). So all 
non-trivial scattering in Eq. (Ilj) stems from J(Q 2 — S 2 ), 
where S" 2 (acting in spin space) do not differentiate be- 
tween orbital angular-momentum channels L directly. 
The quasi-spin operator Q 2 however, do differentiate in 
L, but scatters only in the monopole channel (L = 0) [14]. 
This proves our point; the Slater-Condon interaction, 
at the Laporte-Platt degenerate point F^ /_F( 2 ) = 9/5, 
is free from anisotropic orbital-multipole interactions. 
Here on we will refer to these interactions as simply 
"multipole-interactions" . 

Guided by our findings we propose the following re- 
parametrization of the Slater-Condon interaction 

F<°> = U-U, F^ = 49 (- + I) J, F^ = ? J, (2) 
5 \1 7J 5 

using U, J and 1/7 that control the relative strength 
of the multipole interactions. A cautionary remark 
is in place, the multipole parametrization is arbitrary 
(Ref. [12] uses another equivalent form), and the choice 
of I/7 in Eq. pi) is a matter of taste. 1 However the 
Kanamori limit, without multipole terms, is well defined 
by I/7 = 0. In what follows we set 1/7 = 1/4, which 
corresponds to FW JF^ ss 0.65, in the relevant regime 
for the transition metals [2D] . 

We are now in a position to start our study of the d- 
band model. Much can in fact be learned in the limit of 
strong interactions, where the system turns in to an en- 
semble of isolated atoms with known n-electron ground 
state energies E n . For the ensemble with integer average 
filling h the obvious ground state candidate is the ho- 
mogeneous state with energy E n . But there is also the 
possibility of phase-separated mixtures of atomic states 
with m and ni electrons. In general such a mixture has 
the energy Ei"] n2 = E ni + (E ri2 -E ni )(n-ni)/(n2-n 1 ) 
assuming m < n < n-i- 



1 Our choice is motivated by the simple form Eq. |2J| takes in 
terms of the Racah parameters, A = U — 3J, B = J/7, C = J. 
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FIG. 1: (color online). Ensemble energies as a function 
of J/U for all mixed- valence states (gray lines) relative to 
the atomic d 6 ground state (black line) for 1/7 = 1/4. The 
energy crossings jdi and jdi (dotted lines) in to the d 5 +d 7 
valence-skipping (red line) and d 5 +d 10 split-valence (green 
line) phases are indicated. 
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FIG. 2: (color online). Ensemble phase diagram in the 
(J/U, I/7) and (F (4) /F (2) , F (0) /F (4) ) plane for integer aver- 
age fillings n. The valence-skipping d n_1 +d n+1 phase is only 
present for n = 1,4,6 and 9. The dotted line corresponds to 
l/ 7 = 1/4. 

We have compared all candidate states for every inte- 
ger n and located the ground state crossings as a function 
of J/U and 1/7, see Fig. [ljfor an example. We find that 
(like in Fig. [I]) the valence-skipping d n ~ 1 +d n+1 state is 
the ensemble ground state in the range, jdi < J/U < J42, 
with I/7 dependent bounds jdi = 1/(3 + 8/7) an d 
3d2 = 1/(3 + 2/7), but only for n = 1,4,6 and 9. For 
the other integer n the d n ~ 1 +d n+1 state never becomes 
the ground state. When J/U > jd2, the ensemble has a 
split-valence type of ground state for all n, composed by 
d°+d 5 for n < 5 and d 5 +d w for n > 5 (as in Fig.[l]). The 
final phase diagram is shown in Fig. [2] 

With this background we can understand the connec- 
tion between valence-skipping and effective negative-U. 
The effective Hubbard repulsion U e g is given by [5] 

U e s = E n+ i + E n _i — 2E n = 2{E^ l } l fl+1 — E n ) , (3) 

and L^cff < occurs only for concave series En, 
En + i. In the case of a valence-skipping ensemble ground 
state d n ~ 1 +d n+l we are guaranteed that E < ^ 1 } 1 fi+1 < E n , 
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FIG. 3: (color online). Atomic ground state energy contribu- 
tions for all iV-electron fillings; density-density (squares), spin 
and quasi-spin scattering (triangles), and the multipole ener- 
gies _B' mp ' (circles), up to (irrelevant) linear shifts pN. Note 
that only i?( mp ) is locally concave, and only for JV = 1, 4, 6, 9, 
where the valence-skipping state has a lower multipole con- 
tribution (red lines). 

and Eq. ([3| directly gives U e g < 0. 

But what is now the role of the multipole interactions? 
From the ensemble-phase diagram (Fig. [2]) it is clear that 
the multipole interaction strength I/7 directly controls 
the extent of the valence-skipping phase, and in the limit 
l/7~>0 this phase disappears. We conclude that the 
valence-skipping ground state is realized by the multipole 
interactions. 

To understand why the effect is limited to only certain 
fillings we decompose the atomic ground state energies 
E n in Kanamori and multipole contributions. As seen 
in Fig. [3] the isotropic and mono-pole terms are convex 
(as long as U ~ 3 J > 0). However the multipole energy 
Et p) = E n -(U-3J)N(N-l)/2-J((Q 2 )-(S 2 ))is lo- 
cally concave, but only for the special fillings n = 1, 4, 6, 9 
and can therefore give U e g < for sufficiently large J/7. 
Because of this we will henceforth denote these fillings as 
"multipole-active" . From Fig. [3] it is also clear that the 
valence-skipping d n ~ 1 +d n+1 ground state at multipole- 
active filling n, is stable with respect to doping in the 
whole range n — 1 < n < n + 1 of d n fillings. 

Let us close the discussion of the ensemble by recasting 
the valence-skipping criteria J/U > jdi in terms of 
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As F^'/F^ 2 ' varies weakly within the transition metals, 
fulfillment of Eq. Q is mainly driven by ligand induced 
effective monopole screening of F^ [5]. 

With these insights we leave the strong coupling limit 
and consider the full model with its competition between 
itineracy and local interactions. The ground state is cal- 
culated using the variational Gutzwiller method [2"TH2"3"] , 
previously shown to give phase diagrams in qualitative 
agreement with dynamical mean field theory [21]. Wc 
limit the discussion to translationally invariant param- 
agnetic wave-functions, employing the most general vari- 
ational ansatz with the symmetry of our model. 

Here we report results for d e (particle-hole symmet- 
ric to rf 4 ), whose phase-boundaries are shown in Fig. [4J 
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FIG. 4: (color online). Phase diagram for d 6 , with 1/7 = 
1/4, showing the contours of the local entanglement entropy of 
the metallic state, and the metal-insulator phase-boundaries 
(blue lines), first and second order transitions are indicated 
(solid and dotted lines respectively). 

together with the local entanglement entropy contours 
[25] of the metal. The low J/U region (J/U < j dl ) 
agrees qualitatively with the three-band Kanamori model 
[351 HZ] and w iH not be discussed further. A quantitative 
comparison is left for future works [TS] . 

Our current interest lies in the Hund's-metal [18j and 
valence-skipping regimes (jdi < J/U < jd2)- In gen- 
eral for fixed J/U there is a critical coupling U = U c 
where the metal-insulator transition occurs. But as seen 
in Fig.[4j the U c of the Hund's-metal grows with increased 
J/U, and when J/U — > jdi it diverges (U c — > 00). At 
this point, J/U = jdi, the metallicity prevails for any 
U because the energy-cost for charge fluctuations is zero, 
U e fi(jdi) — 0. When entering the valence-skipping regime 
(jdi < J/U < jd.2), U c becomes finite again as a reen- 
trant valence-skipping d 5 +d7 insulator emerges. Yet, ap- 
proaching the upper boundary J/U —> jd2, U c diverges 
again. Further increasing J/U rapidly reduces U c in favor 
of a split- valence d 5 +d 10 insulator. 

How is then the metal influenced by the change in the 
ensemble ground state from d 6 to d 5 +d 7 ? In terms of lo- 
cal valences the single-particle hopping in the metal gen- 
erates a distribution of adjacent valences. This distribu- 
tion however, is strongly dependent on the intra-atomic 
interaction. 

To investigate this we compute the reduced local many- 
body density matrix p [25] . and calculate the valence 
weights pn as traces of p in every iV-electron subspace. 
The valence-distributions pN for the Hund's and valence- 
skipping metals are shown in FigJEJ at the points marked 
out in the phase diagram (Fig. K In each case pn for 
the corresponding non-interacting metal (U — 0) and in- 
sulator (U — > 00) are shown for comparison. 

The Hund's-metal in Fig. [5^,, has an atomic-like va- 
lence distribution that is substantially narrower com- 
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FIG. 5: (color online). Histograms of valence weights pN as 
a function of filling N for the points marked in Fig. [4] with 
I/7 = 1/4. The correlated metal (red) is shown together 
with the corresponding U=0 non-interacting metal (green), 
and the U — > 00 insulator (blue), for a) the Hund's-metal, 
and b) the valence-skipping d 5 +d 7 metal. 



pared to the non-interacting metal. Most of the weight is 
concentrated in the range N — 5 — 7, with a strong preva- 
lence towards the total average valence N = 6. Turning 
to Fig. and the valence-skipping metal, we find the 
same narrowing down of the distribution, but without 
any certain valence prevalence. Thus, in comparison to 
the Hund's-metal, there is a substantial reduction of the 
average d 6 valence. This type of reduction is the hall- 
mark of the anomalous valence fluctuations driven by the 
effective negative-U in the valence-skipping region. 

We have shown that the local multipole interactions 
drastically reduce the effective Hubbard repulsion U e g in 
the d-band, even making it possible to reach negative- 
U (U c s < 0). Moreover this multipole reduction is only 
obtained for four out of ten possible integer d-band fill- 
ings, namely d 1 , d 4 , d 6 and d 9 . Admittedly we have 
used an over-simplistic model of the d-band. But the 
valence-skipping active fillings is a fundamental property 
of the Coulomb interaction, and apply to the entire class 
of transition metals. 

Experimentally valence-skipping is most clearly ob- 
served when accompanied by charge-order and U c g < 0, 
as in the iron d compounds discussed in the introduction 
[IHZ], and noble-metal d 9 systems such as CsAul3 [28] . 
However multipole-reduced but positive U e g > also 
generate valence-skipping in terms of polarons at elevated 
temperatures T > U e s- This type of thermally induced 
valence-skipping has been used to explain the polaronic 
conduction in d 6 (La,Ca)Co0 3 [29] and d 4 (La,Ca)Mn0 3 
[3"0] . For the d 1 filling pQ some of the candidate transition 
metal complex-oxide compounds are not even thermody- 
namically stable, e.g., La2V20y phase-separate directly 
to LaV0 3 and LaV0 4 (d 1 -> d°+d 2 ) [3T]. 

So returning to the propositions of Anderson [5] and 
Varma [TJ, we conclude that for multipole-active fill- 
ings the electron interaction can drive valence-skipping 
even in absence of lattice relaxation. One such ex- 
ample is Lai/ 2 Cai/ 2 Fe03 [4] that charge-orders to 
3(d 35 ) — > 2(d 3 )+d 5 . While in other cases both multipole- 
interactions and static lattice relaxation combine to give 
[/ c ff < 0. Note that, even thought the rules for multipole- 



active fillings were derived for the degenerate d-band 
model, they remain applicable also in weak crystal- 
fields. However in the limit of strong crystal-fields they 
break down, like in YNi0 3 that show dJ valence-skipping 
charge-order isolated to the crystal-field split e ff -states, 
t%el^t%e° g + t%el 0. 

From the Gutzwiller calculations it is clear that mul- 
tipole interactions also affect the metallic state. The im- 
portance of the metallic valence distribution has been dis- 
cussed in a recent study of SrCo0 3 using the Kanamori 
interaction 32 . Here a followup study also including 
multipole-interactions through the Slater-Condon inter- 
action would be very interesting. 

Although valence-skipping is experimentally most ev- 
ident in charge ordered compounds, negative-U is also 
a potential electron-pairing mechanism for supercon- 
ductivity [TJ]. Therefore it is worth noting that the 
cuprate (<i 9 :Cu 2+ ), ruthanate (d 4 :Ru 4+ ), and iron pnic- 
tide and chalcogenide (<i 6 :Fe 2+ ) superconductors, all 
have multipole-active d-band fillings. 

To conclude we have shown that, in the vicinity of 
the multipole-active fillings d 1 , d 4 , d 6 , and d 9 , the mul- 
tipole part of the Slater-Condon interaction can alone 
drive valence- skipping d n — > d n ~ 1 +d n+1 and negative-U 
in the degenerate d-band. Further more the valence fluc- 
tuations in the valence-skipping metal are drastically dif- 
ferent compared to the atomic-like Hund's metal. None 
of these effects are captured by the Kanamori interaction, 
due to its lack of anisotropic multipole interactions. 
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